February 5, 2008 



Towards absolute scales of radii and masses of open clusters 

A.E. Piskunov^'^'^, E. Schilbach\ N.V. Kharchenko^'^'^, S. Roser\ and R.-D. Scholz^ 



o 
o 

(N 
X) 

O 



> 

(N 
O 
l> 

o 

Oh 

6 



Astronomisches Rechen-Institut, MonchhofstraCe 12-14, D-69120 Heidelberg, Germany 

email: apiskunov@ari.uni-heidelberg.de, elena@ari.uni-heidelberg.de, nkhar@ari.uni-heidelberg.de, roeserfi 

Institute of Astronomy of the Russian Acad. Sci., 48 Pyatnitskaya Str., 109017 Moscow, Russia 

email: piskunov@inasan.rssi.ru 

Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D- 14482 Potsdam, Germany 
email: apiskunov@aip.de, nkharchenko@aip.de, rdscholz@aip.de 
Main Astronomical Observatory, 27 Academica Zabolotnogo Str., 03680 Kiev, Ukraine 
email: nkhar@mao.kiev.ua 



I ari . uni-heidelberg . de 



Received 10 January 2007; accepted , 



Abstract 



Aims. In this paper we derive tidal radii and masses of open clusters in the nearest kiloparsecs around the Sun. 

Methods. For each cluster, the mass is estimated from tidal radii determined from a fitting of three-parametric King's profiles to the 

observed integrated density distribution. Different samples of members are investigated. 

Results. For 236 open clusters, all contained in the catalogue ASCC-2.5, we obtain core and tidal radii, as well as tidal masses. The 
distributions of the core and tidal radii peak at about 1.5 pc and 7-10 pc, respectively. A typical relative error of the core radius lies 
between 15% and 50%, whereas, for the majority of clusters, the tidal radius was determined with a relative accuracy better than 20%. 
Most of the clusters have tidal masses between 50 and 1000 m©, and for about half of the clusters, the masses were obtained with a 
relative error better than 50%. 

Key words. Galaxy: open clusters and associations: general - solar neighbourhood ~ Galaxy: stellar content 



1. Introduction 

As a minimum characteristics to describe a stellar cluster one 
needs to specify the position of a cluster centre and its apparent 
(angular) size. Both parameters serve as a means of identifying 
a cluster on the sky. Also, these parameters are the most com- 
mon, present in numerous catalogues of clusters, and presently 
available for about 1700 galactic open clusters (see e.g. Dias et 
al. I2002i) . In the majority of cases, however, they are derived 
from visual inspection of the area of a cluster on the sky. So, 
the parameters may be strongly biased due to the size of the de- 
tector field, and/or by contamination of field stars, and, hence, 
present a lower limit of the real size of a cluster. As shown by 
KJiarchenko et al. (2005a) these data from literature are normally 
smaller by a factor of two with respect to cluster radii drawn 
from the analysis of a uniform membership based on photomet- 
ric and spatio-kinematic constraints. These latter data are in turn 
subject to various biases (see Schilbach et al. 2006) and need to 
be reduced to a uniform system in order to allow physical in- 
sight into structural properties of the population of galactic open 
clusters. 

Besides the morphological description of a cluster itself, the 
structural parameters carry important information on its basic 
physical properties like mass, and on the surrounding galactic 
tidal field (von Hoerner|T957]). King ( 119621 ) has proposed an em- 
pirical set of cluster parameters as a quantifier of the structure of 
spherical systems, and has shown later (King 1966) that they 
correspond to theoretical density profiles of quasi-equilibrium 
configurations and suit well to stellar clusters with masses span- 
ning from those of open clusters to globular ones. King ( 119621 ) 
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introduced three spatial parameters rc, rt and k, hereafter re- 
ferred to as King's parameters (vc is the so-called core radius, 
rt the tidal radius, and fc is a profile normalization factor), fully 
describing the distribution of the projected density in a cluster 
Since that time, this parameter set is widely used to quantita- 
tively describe the density laws of globular clusters. Here we 
mention a number of studies, where King's parameters were de- 
termined both for galactic globular clusters (e.g. Peterson and 
King [19751 Trager et al. fl^ Lehmann & Scholz IT9971), a nd 
for extragalactic ones (Kontizas 119841 Hill & Zaritski 120061 ) in 
the SMC, or LMC (Elson et al. [T987] ). 

The application of King's parameters might be useful to 
open clusters as well, especially from the point of view of es- 
tablishing a uniform scale of structural parameters and pro- 
viding independent estimates of cluster masses. However, the 
literature on the determination of King's parameters of open 
clusters is much poorer than that of globulars. We mention 
here the following studies based on three-parameter fits of 
selected clusters: King (,1962) . Leonard (1988), Raboud & 
Mermilliod ( I1998al[l998bl ). Among other complicating reasons 
such as an insufficient stellar population, and heavy and irregular 
fore/background, making the study of King's parameters in open 
clusters difficult, we emphasise the difficulty of acquiring data in 
wide-field areas around a cluster. The latter is due to the much 
larger size of a typical open cluster compared to the fields of 
view of detectors currently used in studies of individual clusters. 
Recently, when a number of all-sky catalogues has become avail- 
able, studies exploring the unlimited neighbourhood of clusters 
have been pubHshed (Adams et al. '20011 120021 Bica et al. l2005al 
,2005b 2006). Froebrich et al. (2007) have been searching the 
2MASS survey for new clusters, and provided spatial parame- 



ters for all newly identified cluster candidates derived from the 
fitting of surface density patterns with King profiles. 

Mass is one of the fundamental parameters of star clus- 
ters. There are several independent methods to estimate cluster 
masses. Each of them has its advantages or disadvantages with 
respect to the other ones. But so far there is no method, which 
can be regarded as absolutely satisfactory. The most simple and 
straightforward way is to count cluster members and to sum up 
their masses. Since there is no cluster with a complete census 
of members, one always observes only a subset of cluster stars, 
truncated by the limiting magnitude and by the limited area cov- 
ered by a study, and, therefore, masses from star counts should 
be regarded as lower estimates of real mass of a cluster The 
extrapolation of the mass spectrum to an unseen lower limit of 
stellar masses along with some template of the IMF frequently 
applied in such studies, leads to unjustified and unpredictable 
modifications of the observed mass and should be avoided. The 
farther away the cluster is, the larger is the uncounted fraction 
of faint members, often residing in the cluster periphery. In fact, 
this method could be applied with reasonable degree of safety to 
the nearest clusters observed with deep, wide-area surveys, and 
so providing secure and complete membership. Due to its sim- 
plicity the method is currently widely accepted, and possibly it 
is the only technique which is applied to relatively large sam- 
ples of open clusters (see Danilov & Seleznev |I994| Tadross et 
al. l2002l and Lamers et al. l2005T l. 

The second method is the classical one, namely the appli- 
cation of the virial theorem. It gives the mass of a cluster from 
an estimate of the stellar velocity dispersion, and average inter- 
stellar distances. It does not require the observation and mem- 
bership determination of all cluster stars. The application of the 
method is, however, limited to sufficiently massive stellar sys- 
tems (globulars and dwarf spheroidals) with dispersions of in- 
ternal motions large enough to be measurable. For open clusters 
with typcal dispersions of the order of or less than 1 kms~^ 
present-day accuracies of both proper motions and radial ve- 
locities are fairly rough and are marginally available for a few 
selected clusters only. In spite of this, several attempts have 
been undertaken for clusters with the most accurate proper mo- 
tions (Jones [T9701 [TWll McNamara & Sanders \Wn\ [T9831 
McNamara & Sekiguchi T986 Girard et al. 1989' Leonard & 
Merritt 1989 ), or for clusters with mass determination from ra- 
dial velocities (Eigenbrod et al. l2004l l. 

The third method uses the interpretation of the tidal interac- 
tion of a cluster with the parent galaxy, and requires the knowl- 
edge of the tidal radius of a cluster Considering globular clusters 
which, in general, have elliptical orbits. King (1962) differenti- 
ates between the tidal and the limiting radius of a cluster For 
open clusters revolving at approximately circular orbits one can 
expect the observed tidal radius being approximately equal to the 
limiting one. Though a probable deviation of the cluster shape 
from sphericity may have some impact onto the computed clus- 
ter mass. Nevertheless, this method gives a mass estimate of a 
cluster (Raboud & MermilUod 1 998a|[T998 F) which is indepen- 
dent from the results of the two methods mentioned above. Due 
to the cubic dependence on rt, masses drawn from tidal radii are 
strongly influenced by the uncertainties of rt, however. Taking 
these circumstances into account, one usually reverses the rela- 
tion and calculates tidal radii from counted masses. 

For our studies of open clusters we use the All-Sky Compiled 
Catalogue of 2.5 milHon star^] (ASCC-2.5, Kharchenko ESoTT ). 
including absolute proper motions in the Hipparcos system, B, 

' ftp://cdsarc.u-strasbg.fr/pub/cats/I/280A 
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Figure 1. Apparent density profiles of two open clusters used 
for the determination of tidal radii. The Pleiades are shown 
in the top panel, NGC 129 in the bottom panel. The differ- 
ent colors indicate different samples of stars considered: black 
for l(7-members (group (a), see text), magenta for 2(T-members 
(group (b)), blue for 3(T-members (group (c)), and green for "4(T- 
members" (group (d)). Vertical lines show the radii ri (solid) and 
r2 (dashed). 



V magnitudes in the Johnson photometric system, and supple- 
mented with spectral types and radial velocities if available. 
The ASCC-2.5 is complete down to about V = 11.5 mag. 
Based on the ASCC-2.5 we were able to construct reliable com- 
bined kinematic -photometric membership probabilities of bright 
stars {V < 12) for 520 open clusters (Kharchenko et al. i20D4l 
Paper I), to compute a uniform set of astrophysical parameters 
of clusters, (Kharchenko et al. 12005 al Paper II), as well as to 
identify 130 new clusters (Kharchenko et al. I2005bl Paper III) 
in ASCC-2.5. Currently, we have a sample of 650 open clusters, 
which is complete within a distance of about 1 kpc from the Sun. 
This sample was used to study the population of open clusters in 
the local Galactic disk by jointly analysing the spatial and kine- 
matic distributions of clusters (Piskunov et al. 120061 Paper IV), 
for an analysis of different biases affecting the apparent size of 
open clusters, and the segregation of stars of different mass in 
open clusters (Schilbach et al. i2006i Paper V). 

In this paper we determine King's parameters and tidal 
masses for a large fraction of open clusters from our sample to 
get an independent basis for the construction of a uniform and 
objective scale of spatial parameters and masses. In Sec. |2] we 
briefly discuss our input data, Sec.|3]contains the description of 
the pipeline we apply for the determination of King's parame- 
ters, in Sec.|4]we construct and discuss our sample, in Sec.|5]we 
compare our results with published data on King's parameters 
and with independent estimates of cluster masses. In Sec.|6]we 
summarize the results. 



2. Data 

In this study, we make use of our results on the cluster mem- 
bership (Paper I), and on the parameter determination (Paper II 
and Paper III, respectively) for 650 open clusters. Together with 
other basic parameters like the position of the cluster centre, 
age, distance, and angular size (the apparent radii of the core 
ri and the corona r2), stellar density profiles in the wider neigh- 
bourhood of each cluster are available from our data. According 
to the membership probability PaJE density profiles have been 
constructed for four groups of stars in each cluster: (a) - the most 
probable members {Pap > 61%), so called Icr-members, (b) - 
possible members {Pap > 14%), or 2o'-members, (c) - stars with 
Pap > 1%, or Scr-members, and finally, (d) - all stars in the clus- 
ter area which, for convenience, we call "4cr-members". 

For the construction of the density profiles, we count stars 
in concentric rings around the cluster centre up to 5r2, where 
r2 is the apparent radius of the corona. Due to the relatively 
bright magnitude limit of the ASCC-2.5, the number of clus- 
ter members available for the profile construction is rather low 
(on average, about 45 2cr-members per cluster). In general, the 
number of identified cluster members decreases with increasing 
distance modulus of a cluster. In order to have a statistically rel- 
evant number of stars per concentric ring, one has to chose steps 
of larger linear size for remote clusters. Therefore, we count stars 
in concentric rings of equal angular width (0?05). This homoge- 
neous approach makes the linear profile spacing automatically 
larger for more distant clusters. Since in the following analysis 
we intend to derive 3 unknown King's parameters from the pro- 
file fitting, we need at least four bins in the observed profile, and 
therefore, we select only clusters with r2 ^ 0?2 and with more 
than 10 stars above the background level within r2- Under these 
constraints, the sample contains 290 clusters, though, not for all 
of them, an acceptable solution has been obtained (see Section 
Ell. 

In order to give the reader an idea on the quality of the input 
data for the determination of tidal radii, we show two examples 
of apparent density profiles in Fig[T] the Pleiades, one of the 
"best-quality" clusters in our sample (top panel), and NGC 129, 
a remote cluster, and one of the "low-quality" clusters (bottom 
panel). The corresponding data on the density profiles for all 
650 clusters can be found in the Open Cluster Diagrams Atlas 
(OCDA) available from the CD^ 



3. Determination of tidai radii and ciuster masses 

3. 1. Fitting cluster profiles 

The method we apply is based on the well-known empirical 
model of King (119621) . describing the observed projected den- 
sity profiles /(r) in globular clusters with three parameters Tc, 
rt, and k: 



fir) = fc { [1 + [r/r^r] - [l + [n/r^f] 



(1) 



According to King's definition, is the core radius, rt is the 
tidal radius (approximately equal to the limiting radius in case 



Pap is a newly defined combined membership probability com- 
puted for each star in a cluster area by taking into account only the 
kinematic and photometric criteria, without spatial selection criteria. 
Pap = min{Pfj,, Pph}- Note that this definition of the combined mem- 
bership probability differs from that of Pc used in Paper I. 

' ftp://cdsarc.u-strasbg.fr/pub/cats/l/A+A/438/1163/atlas and 
[ftp:// cdsarc.u-strasbg.fr/pub/cats/J/A+A/440/403/atlas| 



of an open cluster), and fc is a normalization factor which is 
related to the central density of the cluster. This approach has 
also been successfully applied for the determination of the King 
parameters in several nicely populated open clusters (see e.g., 
Raboud & Mermilhod ; 1998al [T998bi ). Nevertheless, the direct 
way of fitting the observed density distribution by the model 
(eq. ([TJ) does not work for the majority of clusters of our sample. 
The main reason is the relatively bright magnitude limit of the 
ASCC-2.5 and, consequently, the low number of cluster mem- 
bers that causes uncertainties in the observed density profiles. 
In order to weaken the influence of poor statistics, we use the 
integrated form of King's formula: 



n(r) 



nrlk |ln[l + (r/r^)^] 



[l + (r/rc)' 



I 1/2 



1 



[1 + in/r,) 



211/2 



(2) 



where n(r) is the number of stars within a circle of radius r. 
Since the spatial boundaries of observed open clusters are not 
clearly defined, and the proportion of field stars projected on 
the cluster area is relatively high, we must take special care of 
choosing the integration limits and the background level if using 
eq. ©. 

Contrary to globular clusters, where the statistical errors of 
empirical density profiles are negligibly low, and the data fix a 
model safely in internal regions of the cluster area, in open clus- 
ters a fit based on the inner area is less reliable and can lead to a 
significant bias in the resulting tidal radius. Therefore, we must 
consider the behaviour of the density profile in exterior regions 
of a cluster and even outside the cluster limits which, a priori, 
we do not know. On the other hand, as one can see from eq. ([T]i, 
the value of f{r) increases at r > rt and tends to a finite limit 
/(r cxd), whereas n{r) goes to infinity for r — ^ cxd in eq. (|2]i. 
This contradicts the physical meaning of n{r) since the num- 
ber of cluster members should be finite, independently of how 
far one expands the counts. Therefore, for a physically correct 
application of eq. Q one should complement it by a boundary 
condition n{r ^ rt) ~ N where N is the number of cluster 
stars. Again, rt and N are unknown. 

In order to overcome the problem, we tried to find a range 
Ar where n{r) is practically flat (see for illustration Fig.|2] left 
panel). This range includes rt, and its length depends on the con- 
centration c — logrt/rc- The range Ar degenerates to a point 
r — rt at c — 0, and increases with increasing c. 

For each cluster, a given membership sample may include 
a number of field stars having, by chance, the proper motions 
and photometry which are compatible with those of real clus- 
ter members. Therefore, before we can try to localize Ar in 
the empirical integrated density profiles, the residual background 
contamination has to be removed. If not, the profiles would in- 
crease steadily with increasing r. This is especially important for 
3(7- and "4(T"-samples which are strongly contaminated by field 
stars, although, it can be essential for Icr and 2fT-samples, too. 

The background correction was done in a uniform way for all 
clusters and membership samples. Assuming that the majority of 
stars at r > r2 are not cluster members, we took the average den- 
sity of a given membership group in a ring r2 < r < 2 r2 to be 
the initial background level in the internal cluster area (r < r2). 
Outside r2 (^2 < r < 5r2), the initial background was first 
set to the observed stellar density in each radial bin. The final, 
smooth background profile over the complete area < ?■ < 5 r2 
was then recomputed as a running average of the initial values 




Figure 2. The general scheme of fitting cluster density profiles. 
Left panel: a theoretical, integrated King's profile for j ^ = 1 
pc, Tt = 20 pc with the fitting range rjn and the tolerance 6 
as indicated. Right panel: the density profiles of the open clus- 
ter Trumpler 3 from the 2tT-sample. The measured profile is 
shown by the short-dashed line. The empirical profile corrected 
for the background (see text) is indicated by crosses, the bars are 
Poisson errors. The vertical (blue) solid and dashed Unes mark 
the empirical parameters ri and r2, respectively. 



with a filter size of 0?55 and a step of 0?05. Finally, the back- 
ground profile was subtracted from the original density distribu- 
tion. An example of an observed profile before and after back- 
ground correction is shown in Fig. |2] right panel. Except for a 
few cases of poor and extended clusters projected on heavy and 
variable background, we obtained reasonable results. 

The general scheme of the profile fitting is shown in 
Fig- Ef left) and is explained in the following: 

Step 1: we start with the determination of an initial value 
r'2 of the fit radius rfn- Assuming the empirical cluster sizes ri 
and r2 to be Vc and rt respectively, we compute r2 from King's 
profile of the concentration class c = log r2 /ri as the distance 
from the cluster centre where the profile does not differ from 
n{r2) by more than an assumed tolerance Q. The tolerance is 
chosen depending on the observation quality of a given cluster, 
and is defined as the average of the Poisson errors of the profile 
outside r2: 9 = {ar{r > ^2)); Cr = ^/n{r). 

Step 2: now we are able to apply a nonlinear fitting routine 
based on the Levenberg-Marquardt optimization method (Press 
et al.. II993I 1 to eq. ^ for r — 0...r2. As initial-guess param- 
eters we chose = ri and — r2, whereas is obtained 
from the solution of eq. (|2|i at r = r2- Empty bins in the dif- 
ferential density distribution were omitted from the integrated 
profile fitting. As a result of successive iterations, we get King's 
parameters together with their rms errors and a x^-value. The 
iterations are stopped when two successive x'^-values do not dif- 
fer by more than lO^'^, and the solution is accepted when the 
number of iterations is less than 100. Then we compute the x^- 
probability function Qix^, v) which can be used as a measure of 
the goodness of fit. For a given degree of freedom v, Q{x^, v) is 
the probability that the difference between the observations and 
the fitted function can be considered random, and their sum of 
the squares is allowed to be greater than x^. According to Press 
et al. (119931 1. the fit can be accepted when Q{x^, v) > 0.1. 



Table 1. Normalized values of King's parameters computed with 
four different membership groups 





Icr 


2cr 


3(7 


4cr 




0.97 ±0.03 


1.00 ±0.02 


1.01 ± 0.02 


1.01 ±0.03 


ft 


1.05 ±0.01 


1.00 ±0.01 


0.98 ±0.01 


0.97 ±0.01 


k 


0.64 ±0.02 


1.02 ±0.03 


1.11 ±0.02 


1.24 ±0.03 


n2 


0.65 ± 0.01 


1.04 ± 0.01 


1.13 ± 0.01 


1.18 ± 0.02 


c 

Or^ 


1.22 ± 0.02 


0.96 ± 0.02 


0.91 ± 0.02 


0.91 ± 0.02 


Srt 


1.32 ±0.04 


0.99 ± 0.03 


0.90 ±0.03 


0.80 ± 0.03 




1.22 ±0.04 


0.94 ± 0.03 


0.90 ±0.02 


0.94 ±0.03 


Q 


1.17 ±0.04 


1.03 ±0.02 


1.00 ±0.02 


0.80 ±0.03 



Although a choice of r fu — f2 provides a fitting of a reason- 
able portion of King's profile to the observations, we consider 
this fit range as the lowest limit of rf it- 

Step 3: In order to check whether a better convergence can 
be achieved, we run the fitting routine (i.e. Step 2) for differ- 
ent rfit ranging from r'2 to 3 r2. If several acceptable solutions 
(i.e. Q(x^j v) > 0.1) were obtained, we selected that one which 
yielded the smallest rms errors in Tc and rt. 

The complete pipeline including the background elimination 
and the profile fitting has been applied to all four membership 
groups. Although the initial guess ^2 of the fit radius is usually 
1.5 r2, the final fit radius rfit turns out to be about 2 r2, and 
the best rt ranges from 1 7-2 to 2 r2. 

However, we must keep in mind that eqs. ([T]i, assume 
spherical symmetry in the spatial distribution of cluster stars, 
whereas a real open cluster is expected to have an elongated form 
with the major axis directed towards the Galactic centre (Wielen 
1974). Thus, depending on the orientation of the line of sight, 
an observer measures a projection rather than the real size of 
a cluster, and the values of rt derived via eqs. ([T]i, give, in 
general, lower limits of the tidal radii. 



3.2. Cluster masses 



According to King (I1962I I. the mass Mc of a cluster at the galac- 
tocentric distance Rq follows the relation 



AA{A- 



B)rl 



G 



(3) 



where A, B are Oort's constants valid for Rq, and G is the grav- 
itational constant (see Standish [l995l l. Since the bulk of our clus- 
ters is located within 2 kpc from the Sun, the linear approxima- 
tion to the velocity field of the disk seems to be reasonable. Thus, 
Oort's constants could be easily expressed by their local values 
Aq = 14.5 ± 0.8 km/s/kpc, Eq = -13.0 ± 1.1 km/s/kpc derived 
in Paper IV for the open cluster subsystem 

A = Aq-Aq SRg 
A-B = A0-B0-2A0 SRg 

where SRq = {Rq ~ Rg,o)I Rofi^ and Rg,o = 8.5 kpc. 

A relative random error of the cluster mass Sm^ = £m^/Mc 
can be derived from eq. (O as 



5^ 



Ma 



9 Si 



2/1 

A-B 



51 



B 



A-B 



Si 



(4) 



where 5^ — Sr 



-'rt/rt, Sa — £a/A, and Sb = Eb/B stand for rel- 
ative random errors of the tidal radius and of Oort's constants, re- 
spectively. The cluster masses and their relative errors are com- 
puted with eqs. (O and (|4|i, and are discussed in Sec. 4.3. 



Table 2. Sample table of King parameters and tidal masses for 236 open clusters. The full table is available in machine readable 
form only. See text for further explanations. 



COCD 


Cluster 


I 


b 


d 


Eb-v 


log* 


ri 


r2 


Us 


is 


^2 


rc 




rt 




k 




log m 








deg 




pc 


mag 


yrs 


deg 










pc 








m 





3 


Blanco 1 


14.17 - 


79.02 


269 


0.01 


8.32 


0.70 


2.90 


4 


2 


53 


1.5 


0.2 


20.0 


3.4 


3.2 


0.4 


3.475 


0.224 


4 


Alessi 20 


117.64 


-3.69 


450 


0.22 


8.22 


0.12 


0.30 


2 


2 


14 


0.5 


0.3 


3.4 


1.3 


10.3 


5.7 


1.152 


0.504 


8 


NGC 129 


120.27 


-2.54 


1625 


0.55 


7.87 


0.11 


0.27 


4 


1 


12 


2.2 


1.1 


8.9 


1.5 


1.0 


0.4 


2.311 


0.226 


11 


Berkeley 4 


122.28 


1.53 


3200 


0.70 


7.08 


0.08 


0.20 


3 


4 


19 


4.3 


2.0 


14.6 


2.7 


0.5 


0.2 


2.815 


0.244 


13 


Alessi 1 


123.26 - 


13.30 


800 


0.10 


8.85 


0.13 


0.45 


4 


1 


19 


1.0 


0.3 


9.1 


2.3 


3.2 


1.3 


2.390 


0.340 






LZ.O.Z.J — 


jU.j I 


zUi 


U.Uj 


cS.j4 


1.70 


1.70 


z 


z 


io 


A A 


1.2 


8.9 


1 1 
1 . 1 


i.Z 


n A 
U.4 


9/111 
Z.41 i 


U. loo 


17 


NGC 457 


126.63 


-4.31 


2429 


0.47 


7.38 


0.12 


0.26 


4 


2 


26 


1.2 


0.6 


13.1 


3.1 


2.6 


1.6 


2.729 


0.310 


22 


NGC 663 


129.46 


-0.94 


1952 


0.78 


7.14 


0.10 


0.30 


4 


1 


24 


3.1 


1.1 


15.3 


3.3 


0.8 


0.2 


2.965 


0.280 


23 


Collinder 463 


127.28 


9.40 


702 


0.30 


8.35 


0.22 


0.72 


4 


3 


77 


3.2 


0.4 


12.4 


1.1 


3.3 


0.3 


2.803 


0.124 


1128 


ASCC 128 


109.93 


-5.96 


900 


0.13 


8.44 


0.15 


0.35 


2 


1 


11 


6.7 


3.5 


7.8 


0.9 


1.9 


1.3 


2.215 


0.156 



Eq. (lUi can be also used to get a "rule of thumb" for the 
prediction of the expected accuracy of cluster masses derived 
from tidal radii. Assuming at first A Aq and B « Bq, one 
obtains Sa ~ 0.06, 5b ~ 0.08. Taking further into account that 
typically > 0.1 (see Sec.|4]l, one finds that the relative error 
of mass determination is dominated strongly by the uncertainties 
of the tidal radius i.e., 

and therefore 

Sm, ~ SSrt- (5) 

4. King's parameters: Results 

4. 1. Construction of the output sample 

For 236 out of 290 clusters in the input list, we obtained at least 
one set of King's parameters by applying the method described 
in Sec. 3.1. Depending on the membership group, each of these 
clusters got from one to four different solutions, and the total 
number of solutions was 708. Per cluster, we have at least one 
set of parameters which are larger than their rms errors, i.e. with 
relative errors of Sr^ — Sr^/rc < 1, — Ert/ft < 1, and 

— £k/k < 1. Since, for the majority of clusters, more than 
one set of King's parameters was obtained, we need a decision 
strategy on the priority of using the results for further analysis. 
From the point of view of membership, the parameters from 1 cr- 
members should be the most reliable, but they are more uncertain 
from the point of view of statistics due to the relative low number 
of stars. The opposite is true for "4o'-members". 

In order to compare the solutions derived from different 
membership groups, to check possible systematics between 
them and to define more or less objective criteria of parame- 
ter selection, we considered a subset of 1 14 clusters. This subset 
includes all clusters having four different solutions. For each of 
these clusters, we computed the mean from the four solutions 
for a given parameter and used this mean as a normalizing fac- 
tor. If a normalized parameter is significantly smaller than one, 
we conclude that a given membership group delivers a signifi- 
cantly smaller parameter than the other groups, and vice versa. 
Table[T|gives the corresponding normalized parameters averaged 
over 114 clusters. 

The most impressive feature of Table [T] is that the core radii 
fc do not depend on the membership group used for computa- 
tion, and tidal radii ft show only a slight systematic dependence 



on the membership groups. On average, the tidal radii obtained 
with the "4(T"-membership samples are smaller only by a factor 
of 1.1 than those with Icr-members. On the contrary, the param- 
eter fc increases towards the "4(T"-membership sample and cor- 
relates strongly with the normalized number of cluster members 
7i2 located within an area of a radius r — r2. This is a logical 
behaviour and follows from the meaning of fc and n-z in eq. (|2]i. 
A relation between the normalized parameter fc and the number 
of the sample i (i — 1, ...,4) used for the solution can be approx- 
imated by 

fc= (0.18 ±0.04) xi + (0. 54 ±0.11). (6) 

As expected,_the relative errors in the determination of the 
parameters Sr^, S^, and Sk are largest for the Icr-solution. The 
l(T-sample contains the most probable cluster members, but their 
number is relatively low compared to the other membership sam- 
ples. This is the reason for relatively large Poisson errors and, 
consequently, for higher rms errors in the fitted parameters. 

The goodness of a fit is given by the Qix^, J^)-probability, 
an output parameter of the fitting pipeline (see Sec. 3.1). The nor- 
malized Q{x^ , t^) -parameters averaged over 1 14 clusters are also 
given in Table [T| As can be seen, Q{x^, v) does not depend on 
the richness of the sample and indicates a more suitable fitting 
with la- and 2cr-samples than with "4cr"-groups. 

Based on the statistics in Table [T] we chose the following 
ranking of the solutions. We give the highest weight to the solu- 
tions with the largest Q{x^, i^) -probability. If, for a given clus- 
ter there are more than one solutions of the same quality (i.e., 
the Q{x^, ;/)-probabilities differ by less than 0.1%), we use an 
additional criterion based on x'^-values supplied by the fitting 
pipeline. Since x^ does depend on the sample size n and on the 
degree of freedom v , a readjusting of the x^-estimate is needed 
when we compare fitting results derived with different mem- 
bership samples. Therefore, we select a solution with a smaller 
value x^n = x'^/i'^ ^ which is, in fact, an average mean 
square deviation of the observed from the fitted profiles in units 
of Poisson errors computed per star In a few cases when even 
the xt „-parameters differ insignificantly (by less than 0.001) for 
two solutions, we give priority to the solution with smaller error 

According to the selection procedure, the solution from the 
l(7-samples gets the best ranking in 94 cases out of 236, from 
2(T-samples - in 72 cases, from Scr-samples - in 39 cases, and 
"4cr"-samples - in 31 cases. 

The data on the structural parameters are compiled in a table 
which is available in machine-readable form only. As an exam- 
ple. Table |2] lists a few entries of the complete data. Column 
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Figures. Distributions of King's radii and their relative errors 
for 236 open clusters. Panels (a) and (b) are for the core radius 
Tc and the tidal radius rt measured in pc. Distributions of rela- 
tive errors in core radius Sr^ and in tidal radius (5,.^ are shown in 
panels (c) and (d), respectively. 



1 gives the cluster number in the COCD catalogue, columns 2 
through 9 are taken from the COCD, while columns 10 through 
20 include the information obtained in this work. For each of 
the 236 clusters we give: name (2), galactic coordinates (3, 4), 
the distance from the Sun in pc (5), the reddening (6), the loga- 
rithm of the cluster age in years (7), empirical angular radii (in 
degrees) of the core ri (8) and of the corona T2 (9). Column (10) 
is the number of the acceptable solutions of King's parameters 
from the four membership groups. Column (11) gives the num- 
ber of the membership sample (la, 2a, 3a, or "4cr") chosen by 
the selection procedure as providing the best solution, whereas 
column (12) is the number of cluster members of this sample 
within r2 after background correction. Columns 13 through 18 
give the corresponding King parameters Tc, rt, A; with their rms 
errors. The parameters are taken as obtained from the selected 
solution. Depending on the applications desired, the reader is 
advised to take into account the empirical relations between the 
selected solution for rt and k and the membership sample used 
to obtain these parameters (cf. Table [T] and eq. Finally, 
columns (19, 20) provide the logarithm of the cluster mass and 
its rms error computed from rt by use of eq. (O and discussed 
in Sec 1431 



4.2. Properties of the sample of King's parameters for 
open clusters 

In Fig. [3a,b) we show the distribution of 236 clusters over the 
derived parameters Tc and rt. The distributions are rather com- 
pact with peaks at about 1 .5 pc for the core radii and at about 
7 - 10 pc for the tidal radii. According to Fig.[3lc,d), the rela- 
tive error Sr^ of the core radius is typically between 15% and 
50%, whereas the tidal radius is more accurate: the majority of 
clusters has the tidal radius determined with a relative error (5^, 



better than 20%. Therefore, for these clusters, we expect to ob- 
tain reasonable estimates of masses from the present data (cf. 
eqs. (|4]i, (|5]l). As a rule, clusters which are more distant and/or 
less populated got relatively large errors in rt. In these cases 
the method reaches its limitations. Probably, if applied to deeper 
photometric data, the method may provide acceptable results for 
a large portion of these clusters. Nevertheless, there is a number 
of clusters with "irregular" density profiles, and, fitting them by 
a model proposed for spherical systems in equilibrium, does not 
have much prospect of success. 

In order to give the reader an idea of typical profiles we show, 
in Fig.m a set of different cases selected from the final sample. 
Compared to the empirical cluster radius r2, the tidal radius rt 
ranges between 1 r2 and 2 r2. For clusters with relatively accu- 
rate radii (Sr^ < 0.33 and (5,.^ < 0.33), the averaged relation 
is 

rt = (1.54 ±0.02) X r2. 

Therefore, we conclude that the empirical parameter r2 scales 
the tidal radius reasonably well. 

4.3. Cluster masses 

Using eq. (O and the tidal radius rt, we estimated masses for 
each of the 236 open clusters. The results are shown in Fig. |5] 
where the distribution of clusters over mass is given in panel 
a), whereas panel b) shows the distribution of relative errors in 
mass. Most of the clusters in our sample have masses in a range 
log Mc/rriQ ~ 1.6 — 2.8, though a few clusters have masses 
as small as IOtoq. Three objects got masses of about IO^toq. 
They are the associations Nor OB5, Sco OB4 and Sco 0B5. For 
139 clusters, the masses were obtained with a relative accuracy 
better than 60%. Their distribution shows the same features as 
the complete sample. 

We note that the masses based on rt from eq. (|2]| do not take 
into account a possible flattening of clusters which arises due to 
the tidal coupling with the Milky Way. Because of this, a stel- 
lar cluster has a shape of a three-axial ellipsoid with the major 
axis oriented in the direction of the Galactic centre. In general, 
we obtain a projection of the tidal radius on the celestial sphere 
from eq. (|2]l, and the relation between the tidal radius and this 
projection depends on the mutual position of the Sun and a given 
cluster Comparing masses of two clusters with different location 
in the Galactic disk, the corresponding effect must, therefore, be 
taken into account. 



5. Comparison with other determinations 

5. 1. King's radii from a three parameter fit 

In contrast to massive spherical systems, there are only a few 
results reported on the determination of structural parameters 
of open clusters via direct parameter fits of King's profiles to 
the observed density distributions. Some of them consider re- 
mote clusters (e.g.. King 1962, Leonard 1988), or newly de- 
tected cluster candidates (Froebrich et al. 2007) which are ab- 
sent in our cluster sample. Others are based on a two parameter 
fit (Keenan[l923l Bica|2006l Bonatto et al. l2005bl) . and, there- 
fore, these results cannot be compared directly with ours. We 
found only seven papers where a three parameter fit was applied 
to open clusters. In these papers only four clusters are in com- 
mon with our sample. These are two nearby clusters, the Pleiades 
and Praesepe, NGC 2168 (M35) at 830 pc from the Sun, and the 
relatively distant cluster NGC 2477, at 1.2 kpc. 
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Figure 4. Examples of radial profiles from our final sample ordered by the priority of the solution. The uppermost row represents the 
Icr-sample. The second row is for 2o--solutions, the third and forth rows illustrate ia and Aa samples. The crosses show empirical 
data corrected for the background, the bars are Poisson errors. The vertical lines indicate the derived radii. The solid lines show r2 
and Tc, the broken lines r2 and rt- The empirical parameters are shown with blue color, while the fitted data are shown with red 
color. Horizontal bars indicate the value of the rms error of the parameter. The solid curve is the fitted King profile, the broken 
curve (yellow) is the total profile not corrected for background. 



In Table[3]we compare the results obtained in this paper with 
those found in the literature for the four clusters. Column 1 is 
the cluster identification. In columns 2 through 6 we provide 
our results: in column 2 we give the number of members within 
the empirical cluster radius r2, column 3 is for Vfn i.e. a radius 
found as best suited for the profile fitting (cf. Sec. 3.1). The re- 
sults of the fitting. King's radii Tc and rt are listed in columns 4 
and 5 together with their rms errors, whereas the corresponding 



tidal mass is given in column 6. Columns 7 through 10 show the 
data from literature: columns 7 and 8 give the number of stars 
used to construct the density profiles and the radius of the area 
considered, respectively. The King radii are listed in columns 9 
and 10. Note that we are not able to keep a uniform format for 
these data due to the different presentation of the results in dif- 
ferent papers. 



Table 3. Comparison of our results with literature data on King's parameters for the clusters in common 
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Present paper 








Literature 
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note 


ref 
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mo 




deg 
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Pleiades 


219 


14 


1.4 ±0.1 


20.5 ± 1.3 


3107 ±638 


270 


6 


1.4 ±0.5 


16 ±7 
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1067 
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0.9...2.9 


13.1 
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1200 


10 


2.1. ..2.8 


12.4... 14.5 


b 


3 


Praesepe 


154 


9 


0.8 ±0.1 


17.1 ± 1.2 


1806 ± 428 


185 


4 


1.0 ±0.5 


11.1 ±4.9 




4 














1000 


3.8 


3.5 


16 


c 
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NGC 2168 


186 


1.95 


1.8 ±0.2 


17.9 ±2.3 


724 ± 677 




0.47 


1.3...4.3 


8.35. ..oo 




6 


NGC 2477 


10 


0.55 


1.9 ±1.0 


6.5 ± 1.2 


89 ±51 




0.25 


1.8 


8.1 


d 


7 



Notes: (a) 

the distribution of low-mass stars (m < 1 raQ); the original results are given in angular units and are transformed to linear sizes by us; (c) - rc, rt 
are obtained from the distribution of low-mass stars (m < 1 m©); (d) - the area is a square of 0.5 x 0.5 sq.deg. 

References: 1- Raboud & Mermilliod tl998ai, 2- Pinfield et al. ( 1998), 3- Adams et al. ( 12001b . 4- Raboud & Mermilhod a998bb . 5- Adams et 
al. ( i2002t . 6- Leonard & Merritt ( |1989; . 7- Eigenbrod et al. ( i2004< . 
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Figures. Distributions of cluster masses (panel (a)) and 
their relative rms errors(panel (b)). Clusters with accurate 
masses((5M < 60%) are shown as the dark histogram in panel 
(a). 



We found three different estimates of the tidal radius pub- 
lished for the Pleiades. A direct comparison with our re- 
sults can be done only with radii obtained by Raboud and 
Mermilliod (I1998al l who consistently applied the 3-parameter 
fitting technique both to various sub-samples of the Pleiades stel- 
lar population and to the total membership sample. From Table 
[3] we conclude that their findings of 7^ and rt coincide well with 
our results. A slightly different method was applied by Pinfield 
et al. dl 9981 1 to derive a tidal radius for the Pleiades. They used 
differential density profiles, which were constructed for cluster 
members falling in different mass ranges, and tidal radii com- 
puted for each of these sub-samples. The results were used to 
derive the partial and total masses of the cluster The tidal ra- 
dius for the cluster as a whole was then computed from the clus- 
ter mass via eq. (O. A similar approach was used by Adams 
et al. (120011) who applied a 3-parameter fit to a sample of low 
mass (m < 1 m©) members of the Pleiades. Again, the tidal 
radius Vt was derived from stellar mass counts in the cluster 
by eq. (O. The approach by Pinfield et al. ( 119981 ) and Adams 
et al. ( 120011 ) provides a smaller value for rt and larger value 
for Tc (for Pleiades members less massive than 1.2 m©) than a 
direct 3-parameter fitting of the density distribution applied to 
the full sample of cluster members. Nevertheless, taking into ac- 
count that the three studies and ours are based on observations 
of different spatial and magnitude coverage, and on an indepen- 



dent membership evaluation, the agreement between the results 
is quite acceptable. 

We arrive at similar conclusions in the case of Praesepe. 
Our estimate for rt is compatible with the result by Raboud 
& Mermilliod (I1998bl) within the rms errors, it coincides well 
with the finding of Adams et al. (120021) who, as in the case of 
the Pleiades, fitted the King profile to low-mass (m — 0.1 — 
1 toq) cluster members. With respect to rc, we achieved good 
agreement with the estimate by Raboud & Mermilliod (1998b| 
whereas rc from Adams et al. (2002) is significantly higher. Such 
a systematic difference can possibly be explained by the consid- 
erably deeper survey used by Adams et al. (2002) in studies of 
the Pleiades and the Praesepe. Since their empirical profiles are 
based on USNO POSS I E and POSS II F plate scans, and on 
the 2MASS catalogue, the input data are dominated by lower 
mass stars which, due to the energy equipartition process, show 
a wider distribution than the more massive stars. 

For the last two clusters in Table [3] there are only two pa- 
pers reporting the determinations of King's parameters. Leonard 
& Merritt ( 1989) studied the central area of the NGC 2168 clus- 
ter Therefore, they were able only to set probable limits for the 
cluster radii. Our estimates of rc and rt fit these limits well. For 
the relatively distant and rich cluster NGC 2477, the King radii 
were published by Eigenbrod et al. ( 120041 ). The authors defined 
the membership sample on the basis of radial velocities of nu- 
merous red giants and constructed density profiles for groups of 
stars of various masses. The core and tidal radii of the cluster 
were determined from the comparison of the structure parame- 
ters of single groups. Also in this case, their results are in good 
agreement with our estimates of rc and rt, though NGC 2477 be- 
longs to the poorer clusters in our sample due to its large distance 
from the Sun. With 10 members within r2, NGC 2477 satisfies 
marginally our constraints. Nevertheless, the fitted parameters 
coincide well with the corresponding estimates obtained with the 
much deeper survey {V < 17) by Eigenbrod et al. ( 120041) . 

5.2. Various scales of cluster masses 

Although the adopted values of Oort's constants and of a dis- 
tance to a cluster can slightly influence its mass estimate, the 
tidal radius is the major source of uncertainty in the mass de- 
termination via eq. (O, simply due to the cubic power rela- 
tion between tidal radius and cluster mass. Analysing different 
groups of the Pleiades members, Raboud & Mermilliod ( 1998al) 
concluded that the tidal mass of the Pleiades is about 14OOm0 
with a 1(7 confidence interval of [530, 2900] niQ. With Oort's 



constants A ^ 15 km/s/kpc, B = —12 km/s/kpc and a clus- 
ter distance of 125 pc adopted by the authors, a strict ap- 
plication of eq. (O would provide a tidal mass of 162Om0 
([330, 4600] m©) for the Pleiades. From a similar analysis 
for Praesepe, Raboud & Mermilliod (1998b) derived a tidal 
mass of 440 7710 ([157, 987] trq) (a consequent use of eq. (O 
would give 520777© ([90, 1570] 7710)). We conclude that the 
disagreement between our estimates and those of Raboud & 
Mermilliod ( |1998al |1998bl l for the Pleiades and Praesepe are 
mainly caused by statistical uncertainties in the determination 
of the tidal radii (cf. Table [Sj. 

In the case of NGC 2477, the situation is less clear. For this 
cluster, Eigenbrod et al. (120041) estimated the tidal and virial 
masses to be Alud ~ 5400 toq and Alyir = 5300 777©, re- 
spectively. This result differs considerably from our estimate. 
Moreover, it is in contradiction to the value of their tidal radius 
of Tt = 8.1 pc which according to eq. ([5]), should give a tidal 
mass of Mtid ~ 170 to©. Therefore, the coincidence achieved 
by Eigenbrod et al. (I2004I I between the tidal and virial masses 
should be considered with caution. We note that NGC 2477 is 
still badly studied, and its membership is poorly established. 
Therefore, a contamination of the stellar sample from Eigenbrod 
et al. (j2t)04) is well probable. This could be one reason for an 
overestimation of the velocity dispersion and, consequently, of 
the virial mass of the cluster. The presently large uncertainties 
in kinematical data for determing virial masses can also cause 
these discrepancies, (see Appendix lAl for more detail). On the 
other hand, this large disagreement can be partly explained by 
an underestimation of the tidal radius: if this relatively distant 
cluster is subject to strong mass segregation, low mass stars on 
the cluster edges can be beyond the magnitude limit even in a 
deep survey. 

Since literature data on tidal masses of open clusters are 
rather scarce, we looked for recent publications on cluster 
masses estimated with other methods. We omit here a discus- 
sion on the determination of virial masses of open clusters and 
refer the reader to Appendix lAl where this method is discussed 
in more detail. In order to compare our results on cluster masses, 
we consider only those publications where cluster masses are ob- 
tained for a relevant number of open clusters rather than for a sin- 
gle cluster. Under these constraints, we found three publications 
on mass determination for galactic open clusters, Danilov & 
Seleznev (fT994b . Tadross et al. (2002), and Lamers et al. (I2005T) . 
In Fig|6]we compare the different results on mass estimates. 

Danilov & Seleznev ( 119941 ) derived masses for 103 compact, 
distant (>1 kpc) clusters from star counts down to S = 16 from 
homogeneous wide-field observations with the 50-cm Schmidt 
camera of the Ural university. For each cluster, the authors esti- 
mated the average mass of a star observed in a cluster and then 
computed the total visible cluster mass. The average mass is 
found either from star counts, or from an extrapolation of the 
Salpeter IMF down to the magnitude limit B = 16. On one 
hand, their cluster masses should be underestimated due to their 
magnitude limited survey, and the bias should increase with in- 
creasing distance modulus of a cluster. On the other hand, with- 
out membership information, the masses could be overestimated 
for clusters located at relatively low distances from the Sun. In a 
certain respect, these biases may partly compensate each other. 

Based on UBV-CCD observations compiled from the litera- 
ture, Tadross et al. ( 120021 ) redetermined ages and distances for 
160 open clusters, and derived cluster masses from counts of 
photometrically selected cluster members. Since they used ob- 
servations taken with different telescopes, i.e., for different clus- 
ters one expects different limiting magnitudes, it is rather diffi- 
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Figure 6. Comparison of tidal masses Mc of open clusters de- 
rived in the present study with literature values on counted- and 
MF-scaled masses as a function of the distance modulus. From 
top to bottom we show ratios of: counted masses of Danilov 
& Seleznev (1994) {Md to tidal masses); counted masses of 
Tadross et al. (120021 ) (Mt to tidal masses); and the IMF-scaled 
masses of Lamers et al. (120051 ) (Ml to tidal masses). The last 
panels show intercomparisons of literature data. 



cult to estimate possible biases. In any case, a large portion of 
clusters should get underestimated masses due to the relatively 
small area of the sky usually covered by CCD observations in 
the past. 

Lamers et al. (2005 ) used data from the COCD for the deter- 
mination of cluster masses. For each cluster, the authors normal- 
ized the Salpeter IMF in the mass range of stars present in COCD 
(V < 11.5), and the normalized Salpeter IMF was extrapolated 
from large masses down to 777 = 0.15 777q. For distant clusters 
with V — My > 8, the extrapolation was done over a very large 
range, from masses larger than « 1.5 7r7o (or My < 3.5) to 
masses of 0.15 7(7,© (or My ~ 13), where the IMF is still not 
very well known. If the IMF at low masses were flatter than 



the Salpeter IMF (see e.g. Kroupa et al. ll993l l. the approach by 
Lamers et al. (I2005I I would give overestimated cluster masses. 
Moreover, in a segregated cluster one should expect different 
forms of the mass function in the central area and at the edges. 
Integration of the Salpeter IMF over the complete cluster area 
would, also, result in overestimating the cluster mass. A com- 
parison of our cluster masses with those of Lamers et al. (I2005I I 
is especially interesting: since both papers use the same obser- 
vational basis, we can estimate uncertainties caused by the dif- 
ferent approaches. 

According to our determination of cluster masses based on 
tidal radii, we expect possible biases for relatively distant clus- 
ters where we observe only a tip of bright stars. Especially in 
segregated clusters, due to energy equipartition the brightest 
stars are more concentrated to the cluster centre, and do not re- 
produce the correct tidal radius. Nevertheless, for distant clusters 
we can see a distance-dependent effect only in the upper panel 
of Fig|6] where we compare our masses with those of Danilov 
and Seleznev ( 1994). Since the A/^/Afc-relations decrease with 
distance modulus, we suppose that our method manages bi- 
ases better than the approach by Danilov and Seleznev (11 9941 ) . 
Moreover, the dependence of the Mu /M^ -relations on distance 
modulus may be explained alone by the biases in the mass deter- 
mination by Danilov and Seleznev (1994) described above, i.e. 
overestimated masses for clusters at low distances and underesti- 
mated masses for distant clusters. This interpretation seems to be 
plausible, but we cannot exclude completely distance-dependent 
biases in our determinations. 

Furthermore, due to a possibly elongated form of open clus- 
ters (cf. Sec. O, our estimates provide a lower limit for cluster 
masses. This can be one of the reasons for systematic differences 
to cluster masses from Lamers et al. (2005 ). Except for nearby 
clusters with V — My < 8, their masses are, on average, larger 
by a factor of 10. On the other hand, masses from Lamers et 
al. (2005J are systematically larger than masses from Danilov 
and Seleznev CI 994) and Tadross et al. (2002^, too. Therefore, 
we can not exclude that the approach by Lamers et al. (120051 ) is, 
at least partly, responsible for these differences. 

6. Conclusions 

In an ideal case, i.e. if, for a given cluster, the membership is de- 
termined with certainty and completeness, the cluster mass could 
simply be derived from counting masses of individual members. 
In the future, with deeper surveys and an increasing accuracy of 
kinematic and photometric data, this primary method will pro- 
vide sufficiently accurate and uniform mass estimates for a sig- 
nificant number of open clusters in the Galaxy. At present, how- 
ever, there is no way to measure the masses of open clusters 
directly. The methods currently applied require a number of as- 
sumptions, and depending on the method and assumptions used, 
the results can differ by a factor of 100 for individual clusters 
(cf. Fig|6l). Therefore, the determination of cluster masses is still 
a very challenging task. 

Our aim was to estimate masses for a larger number of clus- 
ters by applying a uniform and possibly objective method, and 
to obtain an independent basis for statistical studies of the dis- 
tribution of cluster masses in the Galaxy. In our work we could 
benefit from the homogeneous set of cluster parameters derived 
for 650 open clusters with good membership based on the as- 
trometric and photometric data of the ASCC-2.5. The estima- 
tion of cluster masses was done via tidal radii determined from 
a three-parameter fit of King's profiles to the observed density 
distribution (King [1962]). This method is weakly dependent on 



assumptions and can be applied to all spherical systems in equi- 
librium with well-defined density profiles. Since these require- 
ments are not always met in the case of open clusters, we could 
obtain solutions only for 236 clusters, i.e. for less than half of 
the clusters in our sample. However, this number is consider- 
ably larger compared to the small number of clusters with tidal 
masses determined before. 

The main difficulties in the practical application of King's 
model to open clusters arise from the relatively poor stellar pop- 
ulation (compared to globular clusters) and from the higher de- 
gree of contamination by field stars in the Galactic disk. Using an 
all-sky survey, we could rely on the completeness of data in the 
selected sampling areas, down to the limiting magnitude of the 
ASCC-2.5. Further, since we were free in selecting the size of 
the sampling areas, we were able to optimize the boundary con- 
dition for each cluster as well as possible. As a result, we could 
partly decrease the influence of the above mentioned problems 
in applying King's method and improve the solutions by taking 
into account the outermost regions of the clusters and by exclud- 
ing residual field stars from the solution. 

Together with a realistic membership based on both kine- 
matic and photometric constraints, a good profile fitting could 
be achieved even for clusters with a relatively low number of 
members. The highest quality of the fitting (goodness-of-fit) 
was achieved with the best-determined membership sample (so- 
called "lcr"-members) and, hence, a low contamination by field 
stars. However, it turned out that the membership criterium alone 
did not have very strong impact onto the values of the fitted pa- 
rameters Tc and rt themselves. In fact, for compact and relatively 
distant clusters, we sometimes found the best results without a 
preliminary membership selection (so called "4(T"-members). In 
conclusion, this paper could be seen as a justification for a sim- 
ple application of King's method to observed brightness profiles 
of compact open clusters no matter if membership is determined 
or not, provided that the observed density profiles are properly 
corrected for the background. 
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Appendix A: Viriai masses of open clusters: 
current status 

All our attempts failed to compute reasonable cluster masses 
from the dispersion of proper motions and/or radial velocities 
taken from the ASCC-2.5 catalogue for cluster members. The 
main reason is the accuracy of kinematical data, which is still 
too low in current all-sky surveys. 

Up to now, the best published data on velocity dispersions 
were obtained for a few open clusters from proper motions ob- 
tained from long-term observations with the Yerkes 40-inch re- 
fractor (F = 19.3 m, a scale of 10.7 arcsec/mm, a typical epoch 
difference of more than 55 years, and a typical rms-error of the 
proper motions of a few 0. 1 mas/y). The clusters are the Pleiades 
(Jones 197 0), Pra esepe (Jones \Wri\i . NGC 6705 (McNamai-a 
& Sanders [1977] ). NGC 6494 (McNamara & Sanders 1983|, 
NGC 2168 ( McNamara & Sekiguchi 1986j), NGC 2682 (Girard 
et al. 119891 ). In these papers, the internal proper motion dis- 
persions are corrected for different biases, and viriai as well 
as counted masses are determined for four clusters i.e., for 
the Pleiades, Praesepe, NGC 6494, NGC 6705. The results are 
shown in Fig lA.ll 

Independent of the methods of mass determination, the 
masses in Fig. lA.ll show a correlation with cluster distance. 




Figure A.l. Masses and internal velocity dispersions for the 
open clusters: Pleiades, Praesepe, NGC 6494, NGC 2168, 
NGC 2682, and NGC 6705 (the clusters are ordered with dis- 
tance). Upper panel: internal velocity dispersion versus distance. 
Dots are for proper motion dispersions, triangles for tangential 
velocity dispersions. Bars are rms-errors of the dispersions: bars 
with hats are for proper motions and bars without hats for tan- 
gential velocities (for the Pleiades and Praesepe, the symbols 
for the proper motion dispersions are shifted to the left for a 
better visibility). The line shows a linear fit to the tangential 
velocity dispersion data. Lower panel: Cluster mass versus dis- 
tance (based on literature data, see text for more detail). Crosses 
mark virial masses obtained from proper motion dispersions, and 
squares show counted masses. 



Since the number of clusters is very low, the correlation does 
not need to be real, but can appear by chance, due to the 
small sample. In order to understand the effect, we transformed 
the proper motion dispersions cr^ published for 6 clusters to 
one-dimensional tangential velocity dispersions (t„ via (t„ = 
4.74 d (7^, where d is the distance of a cluster adopted in the orig- 
inal papers. In the upper panel of Fig. lA.ll we show cr^ and (t^, 
as functions of distance d. Whereas (t^ is independent from the 
cluster distance, the tangential velocity dispersion indicates a 
strong correlation which is well described by a first-order poly- 
nomial. Therefore, we suppose a non negligible random compo- 
nent in CTu resulting rather from residual rms errors in proper 
motions than from the internal velocity dispersion. This is not 
surprising, because, for all clusters at distances larger than the 
Pleiades and Praesepe, the dispersions are of the order of the 
rms errors of the proper motion measurements. Assuming that 
the internal velocity dispersion were similar for all these clus- 
ters, we obtained ay = 0.31 km/s by extrapolation of the linear 
regression polynomial to d = 0. In other words, we need tangen- 
tial velocities determined with an accuracy better than 0.3 km/s 



(or proper motions with an accuracy better than 0.06-0.07 mas/y 
for a cluster at d = 1 kpc) for a more or less reliable estimate of 
its virial mass. Even the Hipparcos proper motions with typical 
rms errors of 1 mas/yr do not meet this requirement. 

With respect to radial velocities, one may suppose that these 
give a better basis for the determination of the internal veloc- 
ity dispersion since their accuracy does not depend on distance. 
Recently, radial velocities have been measured in several open 
clusters for a sufficient number of cluster members (see e.g., 
Eigenbrod et al. 2004 for NGC 2477, or Furesz et al. l2056l 
for NGC 2264). The derived velocity dispersions are, however, 
somewhat too large, 0.93 km/s in NGC 2477, and 3.5 km/s in 
NGC 2264. This may be a consequence of biases which still af- 
fect the data and which are very difficult to take into account. 
Among them may be the contamination by field stars and by un- 
resolved binaries, or motions within stellar atmospheres. 
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